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ABSTRACT 

We examine the nature and role of mass loss via an equatorial decretion disk in massive stars with near-critical rotation 
induced by evolution of the stellar interior. In contrast to the usual stellar wind mass loss set by exterior driving from 
the stellar luminosity, such decretion-disk mass loss stems from the angular momentum loss needed to keep the star near 
and below critical rotation, given the interior evolution and decline in the star's moment of inertia. Because the specific 
angular momentum in a Keplerian disk increases with the square root of the radius, the decretion mass loss associated 
with a required level of angular momentum loss depends crucially on the outer radius for viscous coupling of the disk, 
and can be significantly less than the spherical, wind-like mass loss commonly assumed in evolutionary calculations. 
We discuss the physical processes that affect the outer disk radius, including thermal disk outflow, and ablation of the 
disk material via a line-driven wind induced by the star's radiation. We present parameterized scaling laws for taking 
account of decretion-disk mass loss in stellar evolution codes, including how these are affected by metallicity, or by 
presence within a close binary and/or a dense cluster. Effects similar to those discussed here should also be present in 
accretion disks during star formation, and may play an important role in shaping the distribution of rotation speeds 
on the ZAMS. 
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1. Introduction 

. Classical models of stellar evolution focus on the domi- 
' nant role of various stages of nuclear burning in the stel- 
I lar core. But in recent years it has become clear that 

■ stellar evolution, particularly for more massive stars, can 
' also be profoundly influenced by the loss of mass and an- 
, Kular momentum from the stellar envelope and surface 

■ (jMaeder fc Mevne^l2008[ ). In cool, low-mass stars like the 
[ sun, mass loss through thermal expansion of a coronal wind 
. occurs at too-low a rate to have a direct effect on its mass 

■ evolution; nonetheless the moment arm provided by the 
] coronal magnetic field means the associated wind angular 
. momentum loss can substantially spin down the star's ro- 
tation as it ages through its multi-Gyr life on the main se- 

, quence. Except in close binary systems, the rotation speeds 
• of cool, low-mass stars are thus found to decline with age, 
' from up to ~100 km s""'^ near the ZAMS to a few km s""'^ 
for middle-age stars like the sun. 

By contrast, in hotter, more massive stars the role 
and nature of mass and angular momentum loss can be 
much more direct and profound, even over their much 
shorter, multi-Myr lifetimes. While some specific high-mass 
stars appear to have been spun dow n by strongly mag- 
netized stellar winds (e.g HP 19161 2, iDonati et al.l[2006l 
or HD 37776. iMikulasek et al]|2008[) . most massive stars 
are comparitively rapid rotators, with typical speeds more 
than 100 km s~^, and in many stars, e.g. the Be stars, 
even approaching the critical rotation rate, at which the 
centrifugal acceleration at the equatorial surface balances 



Newtonian gra vity ()Howarthl [200l iTownsend et al.ll2004 
lHowarthl[2007h . 

Indeed, models of the MS evolution of rotating massive 
stars show that, at the surface, the velocity approaches the 
critical velocity. This results from the transport of angu- 
lar momentum from the contracting, faster rotating inner 
convectivc core to the expanding, slowed down radiative 
envelope (jMevnet et al.l (2006) . In stars with moderately 
rapid initial rotation, and with only moderate angular mo- 
mentum loss from a stellar wind, this spinup from inter- 
nal evolution can e ven bring the star to critical rotation 
(jMevnet et al.l[2007l) . Since any further increase in rotation 
rate is not dynamically allowed, the further contraction of 
the interior must then be balanced by a net loss of angu- 
lar momentum through an induced mass loss. In previous 
evolutionary models, the required level of mass loss has typ- 
ically been estimated by assuming its removal occurs from 
spherical shells at the stellar surface (jMevnet et al.|[2006ll . 

This paper examines the physically more plausible 
scenario that such mass lo ss occurs throug h an equato- 
rial, viscous decretion disk (jLee et al.l 119911) . Such decre- 
tion disk models have been extensively applied to analyz- 
ing the rapidly (and possibly near-critically) rotating Be 
stars, which show characteristic Balmer emission thought 
to originate in geometrically thin, warm, gaseous disks in 
Keplerian orbit near the equatorial plane of the paren t star 
(jPorter fc Riviniud l2003l ICarciofi fc Biork manll2008[ ). But 
until now there hasn't been much consideration of the role 
such viscous decretion disks might play in the rotational 
and mass loss evolution of massive stars in general. 
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As detailed below, a key point of the analysis here is 
that, per unit mass, the angular momentum loss from such 
a decretion disk can greatly exceed that from a stellar wind 
outflow. Whereas the angular momentum loss of a nonmag- 
netized wind is fixed around the transonic point very near 
the stellar surface, the viscous coupling in a decretion disk 
can transport angular momentum outward to some outer 
disk radius i^out, where the specific angular momentum is 
a factor \/ i?out / ^oq higher than at the equatorial surface. 
For disks with an extended outer radius -Rout ^ ^oq, the 
angular momentum loss required by the interior evolution 
can then be achieved with a much lower net mass loss than 
in the wind-like, spherical ejection assumed in previous evo- 
lution models. 

For a given angular momentum shedding mandated by 
interior evolution, quantifying the associated disk mass loss 
thus requires determining the disk outer radius. For ex- 
ample, in binary systems, this would likely be limited by 
tidal interactions with the companion, and so scale with 
the binary separation (jOkazaki et al.|[2002| ). But in single 
stars, the processes limiting this outer radius are less ap- 
parent. Here we explore two specific mechanisms that can 
limit the angular momentum loss and/or outer radius of 
a disk, namely thermal expansion into supersonic flow at 
some outer radius, and radiative ablation of the inner disk 
from the bright central star. For each case, we derive simple 
scaling rules for the required disk mass loss as a function 
of assumed stellar and wind parameters, given the level of 
interior-mandated angular momentum loss. 

The organization for the remainder of this paper is as 
follows: Sect. [5] presents simple analytical relations for how 
the presence of a disk affects the mass loss at the critical 
limit. Sect. [3] develops set of equations governing structure 
and kinematics of the disk, while Sect. U solves these to 
derive simple scaling for how thermal expansion affects the 
outer disk radius and disk mass loss. Sect. [5] discusses the 
effects of inner-disk ablation by a line-driven disk wind in- 
duced from the illumination of an optically thick disk by 
the centeral star, deriving the associated ablated mass loss 
and its effect on the net disk angular momentum and mass 
loss. Sect.Elgives a synthesis of the different cases discussed 
here and offers a specific recipe for incorporating disk mass 
loss rates into stellar evolution codes. Sect. [7] discusses some 
complementary points (e.g. viscous decoupling, tidal effects 
of nearby stars, reduced metallicity, etc.), while Sect.[8]con- 
cludes with a brief summary of the main results obtained 
in this work. 



wind, etc. is negligible, i.e. J = 0, then the star has to spin 
up at a rate given by 



(2) 



n _ I 

However, once the star reaches the critical rotation fre- 
quency Q = f^ciit = 



GM/Rl^ (where M is the steUar 
mass and i?eq is the equatorial radius when the star is ro- 
tating at the critical limit), this spin- up has to end (il = 0), 
requiring instead a shedding of angular momentum given 
by 

j = merit. (3) 

If we assume this occurs purely through mass loss at 
a rate M through a Keplerian decretion disk, the angular 
momentum loss is set by the outer radius i?out of that disk, 
given by 



JK(i?out) = M«K(i?out)i?out = A/r!eriti?e ' ' 




(4) 



cq 



where the Keplerian velocitjQ is vyi{r) — ^ GM/r. Setting 
JK(i?out) equal to the above J required by a moment of 
inertia change /, we find the required mass loss rate is 



M = 



^''oq 



i?, 



cq 



Rn 



(5) 



As i?out gets larger, note that the required mass loss rate 
gets smaller. 

Eq. ^ can be compared with the case where mass de- 
couples in a spherical shell just at the surface of the star, 
i.e., where i?out = ^cq- In that case the required mass loss 
is just 



2RL 



(6) 



So when a Keplerian disk is present, the mass loss is re- 
duced by a factor § \/ Rout / Rcq with respect to the case 
with no disk. If i?out is small, large mass-loss is necessary 
to shed the required amount of angular momentum to keep 
the rotation frequency at its critical value. In the oppo- 
site case, when i?out becomes substantial, only a relatively 
small mass-loss is necessary to shed the required amount of 
angular momentum, implying a significant difference in the 
mass loss evolution. 



2. Basic analytic scaling for disk mass loss 

Let us begin by deriving some simple analytic expressions 
for the effect of the disk viscous coupling on the disk mass- 
loss rate. 

Assuming a star that rotates as a rigid body, the magni- 
tude of stellar angular momentum J is given by the product 
of the stellar moment of inertia / and the rotation angular 
frequency n, J — Ifl. During stellar evolution, the time rate 
of change of angular momentum depends on the changes in 
moment of inertia and rotation frequency. 



J + m. 



(1) 



If, for example, the moment of inertia declines at a rate 
/, and the change of the angular momentum through any 



3. Numerical models 

To obtain a detailed disc structure, we solve stationary hy- 
drodynamic equations in cylindrical coordinates integrated 
over t he height above the equatorial pla ne z (iLightmanI 
\l974, iLee et all Il991l [Okazaki 2001, Jo nes et al.l I2008D . 
Assuming axial symmetry, the corresponding variables, i.e., 
the radial and azimuthal velocities Vr, v^j,, and the inte- 
grated disk density E = J^^pdz, depend only on radius 
r. The equation of continuity in such a case is 



1 d (r'Evr 
r dr 



= 0. 



(7) 



^ When the star is rotating at the critical limit, the critical 
rotational velocity is equal to the Keplerian velocity at i?oq, 



J. Krticka, S. P. Owocki & G. Meynet: Mass and angular momentum loss via decretion disks 



3 




r/R, 



eq 



1000 



10^ 
10 

1 



10 



10 



10"' 



10- 



-2 



10 



-5 





).......•■ 






/ ' 




/ ''' 






p = 




p=0.2 - 




p = 0.4 



10 



10" 



10'' 



10^ 



r/R, 



eq 



10" 



Fig. 1. The dependence of the radial velocity, azimuthal velocity, and the angular momentum loss rate in units of equator 
release angular momentum loss rate JK(-Req) on the radius in a viscous disk. Left: models of isothermal disk {p = 0, 
Tq — ■^Tcs) with different viscosity parameter a. Right: Models with various temperarure profile for fixed a = 0.1 and 
To = ^Tcff. Arrows denote the location of critical points. 



The stationary conservation of the r component of momen- 
tum gives 



dvr 
dr 



r 



1 d(a^S) 
S dr 



3a^ 
2 r ' 



(8) 



where g = —GM/r 



kT/{^mii), with the tem- 



perature assumed to vary as a power-law in radius, T — 
To (Rcq/r)^ , where Tq and p are free parameters, ^ is the 
mean molecular weight (taking /i = 0.62), and mn is the 
mass of a hydrogen atom. In the equation of conservation of 
the (/> component of momentum we introduce the viscosity 
term (jShakura fc Sunvaevlll973h parametrized via a 



r dr 



r^S dr ' 



(9) 



and the conservation of the 9 component of momentum 
gives the hydrostatic equilibrium density distribution 



Po exp 



2lp) ' 



H 



(10) 



Note that the equatorial density po is related to the 
vertically integrated disk density via S = ^/2ttpoH. 
Close to the star, detai l ed energy-balance models 
(|Millar fc Marlboroughl FTool iCarciofi fc BiorkmanI [2008h 
show the disk is nearly isothermal with Tq = ^Toff and 
p = 0. But to account for the radial decline of the tempera- 
ture in the outer regions, we also consider here models with 
power law temperature decline, with p > 0. 

The system of equations Eq. (HI)-® has to be supple- 
mented by appropriate boundary conditions. For obtaining 
the value of Vr at the stellar surface r — R^q we use the 



fact that at the critical point with radius i?crit given by the 
condition 



GM 5 

-I- 



crit 



^ "-crit 



2R. 



crit 



do^ 
dr 



(11) 



we should have that Wr(^crit) = a to ensure the fini teness of 
the d erivatives at this point (Eqs. ([7]), dS]), see also lOkazakil 
I2nnih . Thus we chose Vr at the surface such that at i?crit we 
have Vr = a. Note that the radial disk velocity is supersonic 
above the critical point. The value of the azimuthal veloc- 
ity at the stellar surface is equal to the corresponding 
Keplerian velocity. The system of studied hydrodynamical 
equations is invariant for the change of the scale E' = 7S 
(where 7 is constant). Consequently, the equations do not 
provide any constraint for the mass- loss rate M = 2TTrVr^, 
which in our case is obtained from the angular momentum 
loss required by the evolutionary calculations. This provides 
the remaining boundary condition for the column density 
S. Here we treat the disk mass-loss rate as a free parameter. 

For the numerical solution of the system of equa- 
tions Eq. we approximate the differentiation at se- 
lected radial grid and use the Newton-Raphson method 
(e.g., Krticka 2003). The resulting system of linear equa- 
tions is solved using the numerical pa ckage LAPACK 
(http : // www, cs ■ Colorado . edu/~lapack^ [Anderson et al.l 
1999). 

4. Results of numerical models 

The general disk properties do not significantly depend on 
particular stellar parameters. Nevertheless, to be specific, 
for a detailed modelling we selected the stellar parame- 
ters roughly corresponding to evolved massive first star 
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(|Marigo et al.ll200lL lEkstrom et aI1l2008bD T^s 30 000 K, 

M = 50Mq, R = 30Rq. 

The calculated models for different values of a are given 
in Fig.[TJ Close to the star the integration of the momentum 
equation Eq. ^ using the continuity equation Eq. (O 



Note also that adding cooling can substantially increase 
the critical radius and thus the disk angular momentum 
loss. For example, for p = OA the angular momentum loss 
increases by a factor of 10 compared to the isothermal case 
(see Fig. ID). 



= const. 



(12) 



gives linear dependence of the radial velocity on radius in 
isot hermal disks (f or ti^ <C a), ~ r, consequently E '-^ 
(|Okazakill2001f ). Finally, from the momentum equation 
Eq. ([5]), it follows that close to the star the disk rotates 
as Keplerian one, i.e. ~ r~^/^. As a result, the angular 

momentum loss scales as J ~ rv^ ~ r^/^, in accordance 
with Eq. (|3]) . As the disk accelerates in radial direction, Vr 
becomes comparable with a and the term rv,j, dominates 
in Eq. (|12p , consequently the disk is momentum conserving 
close to the critical point, rv^ = const, (see Fig. [T]). 

In the supersonic region from the momentum equation 
Eq. ^ follows the logarithmic dependence of the radial ve- 
locity on radius, v'^ ^ Inr. Consequently, the second term 
in equation Eq. rises and as a result of this may 
become even negative. However, this behaviour is a conse- 
quence of adopted Shakura-Sunyaev viscosity prescription 
which predicts non-zero torque even for shear-free disks, 
and is likely not applicable in the supersonic region. 

A maximum angular momentum loss due to the disk 
is obtained in the case when the disk has its outer edge 
at the radius where J is maximum (see Fig. [T|). Note that 
this value does not significantly depend on the assumed vis- 
cosity parameter a. An estimate of the maximum angular 
momentum loss can be obtained assuming that it is equal 
to the angular momentum loss at the critical point. From 
the numerical models it follows that the azimuthal veloc- 
ity at the critical point is roughly equal to the half of the 
Keplerian velocity (see Fig. [T|), 



V^{Rcrit) ~ ■^'UK(i?crit) 



(13) 



In this case the critical point condition Eq. (jlip yields an 
estimate of the critical point radius 



Rc 



Re 



f Vk{Rc^ 



10 + 4p V fl(-Rcq) 



(14) 



from which the maximum angular momentum loss via the 
disk follows 



Ja(M) « - 



t^K(-Roq) 
10 + 4p\ a(i?eq) 



MvKiReq)Re 



(15) 



5. Radiative ablation 

As the radiative force may drive large amou nt of mass out 
of the hot stars via line-driven wind (see lOwockil |2004 
iPuls et cd] l2008l for a review) it may also effectively set 
the outer disk radius. The radiative force may in this case 
ablate the material from the disk and sustain a radiatively 
driven outflow (Gayley et al. 1999, 2001). In the following 
we give an estimate of the disk wind mass-loss rate, which 
is derived in Appendix. 

The disk outflow may in our case resemble the radiation 
driven win ds from luminous accre t ion disks jP roga et alj 
199 &. .Feldmeier fc Shlosmanl Il999l iFeldmeier et al. 199S 



Pr oea et al. 19991 ). The outflow in these simulations orig- 



inates from the whole disk surface. Consequently, part of 
the stellar outflow is carried outwards by the disk and part 
by the disk wind and the fraction of material carried out 
by the disk wind increases with radius. The disk will be 
in this case truncated at the radius where the material is 
carried away entirely by the wind. As the viscous coupling 
is likely not maintained in the supersonic wind, only the 
ablation of the material from the regions close to the star 
would decrease the effectiveness of braking. 

Let us roughly de termine the mass-loss rat e of suc h disk 
wind. The classical ICastor. Abbott fc KleinI ()1975L here- 
after CAK) stellar wind mass-loss rate estimate 



McAK 



a 



\ — a c 



(FQ) 



1/q-I 



(16) 



where Q and a are line force parameters (see also iGavlevI 
Il995f ). L is the stellar luminosity, and the Eddington pa- 
rameter F = Kei/ (47rGMc), with beeing the Thomson 
scattering cross-section per unit of mass, can be rewritten 
in the term of mass flux from a unit surface. 



m 



F / K^FQ 



1/q-I 



1 



eg 



(17) 



where F is the driving flux and g is local gravitational accel- 
eration. The radiative energy impinging the unit of surface 
parallel to the direction to the star is from geometrical rea- 
sons proportinal to FR/r, where R is the polar radius, and 
F is the radiative flux at radius r. Assuming that the disk 
is optically thick (see Sect. lA.lj) . and all incident radiation 
is directed upward, we can roughly estimate 



F = ^F 
r 



(18) 



In agreement with Fig. [l] comparing the formula Eq. (jl5p 
with analytical estimate Eq. Q the angular momentum 
loss is roughly given by ^ JK(-Rcrit), i-e., it is one half of the 

angular momentum loss of the Keplerian disk truncated Taking g = GM/r^, the total disk wind mass- loss rate is 

then given by an integral of the mass-loss rate per unit of 
the disk surface rh between the equatorial radius R^q and 
the outer disk radius i?out 



at the critical point radius i?crit- The factor ^ comes from 
the fact that the disk is not rotating as a Keplerian one 
at large radii (see Fig. [1]). Hence, the minimum disk mass 
loss rate required for given moment of inertia decline is by 

a factor of about (wK(-Rcq)/a(^oq))^^^~^ lower than in the 
case without a disk. 



Mdw(i?out) = 2 X 27r 



-Ro 



(19) 
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where factor of 2 in Eq. (TT^ comes from the fact that the 
wind originates from both sides of the disk. Inserting the 
mass flux estimate Eq. dTT]) and ^TE\\ we derive 



Mdw(i?out) = McakPi 
where (using substitution x — r/R) 



Ro 



R 



(20) 



-i/"-^dx 



3/2 



1 



(21) 

Assuming the disk wind is not viscously coupled to the 
disk, the total angular momentum loss rate via the disk 
wind is 



O 



2 X 27r 



mV(f,r^ dr. 



(22) 



As the disk wind originates mainly from the regions close to 
the star (with r/R < 10, see Fig. lA.2|) . where the azimuthal 
velocity is roughly equal to the Keplerian one (see Fig.[T]), 
we can assume ~ VK{r) in Eq. and consequently 
the disk wind angular momentum loss rate 



jdw(i?out) -i?^'K(i?)Pi 



R„ 



R 



McAK 



(23) 



is by a factor of Pi {Rout/ R) larger than the angular mo- 
mentum loss due to the CAK wind launched from equator 
of hypothetical critically rotating spherical star with radius 
R. 

A more detailed calculation (see Appendix |^ gives a 
more complicated form of i-'^(a;out) via Eq. (jA.23p 



PiiXont) 



2-K q32q 



1 



(24) 



which shah be used in Eqs. (gOl), (ESI) instead of Eq. (EH)- 
For an infinite disk (i?out ~^ oo) we derive from Eq. (j24p 
maximum disk wind mass-loss rate 



Mdw(oo) = 2i+i7r-i3-5t-iQ,AfcAK, 
and maximum angular momentum loss rate as 

22 + c7r q3 2q 



(25) 



-'"dw(oo) = 



a 



-aRvK{R)McAK- (26) 



For a typical value of a w 0.6 (|Krtickall2006f) the maximal 
disk wind mass-loss rate is relatively low, just about 1/25 
of the CAK stellar wind mass-loss rate. 



6. Mass loss of the star-disk system at the critical 
limit 

The structure of the decretion disk and the radiatively 
driven wind blowing from its surface depends on the value 
of the angular momentum loss J needed to keep the stellar 
rotation at or below the critical rate and on the magni- 
tude of the radiative force. If the angular momentum loss 
is small, then the disk could be blown away by the radiative 
force already very close to the star. In the opposite case, if 



the angular momentum loss is large, then the mass carried 
away by the disk wind is negligible. 

In the intermediate case the mass and angular momen- 
tum will be carried partly by the disk and partly by the 
disk wind. When the star has to lose angular momentum 
at a rate J to keep at most the critical rotation, the angular 
momentum will be carried by the stellar wind (Jw), by the 
disk wind (Jdw), and by the disk itself (J5), 



J — Jw + Jdw + Ja{Md). 



(27) 



For the calculation of total disk mass-loss rate the following 
procedure could be used. 

For a given stellar and line-force parameters {Q and a) 
the maximum disk wind angular momentum loss Jdw (00) 
corresponding to infinite disk i?out ^ 00 can be calculated 
using Eqs. (j26l) . If the net angular momentum loss that 
should be carried away by the disk outflow J — is lower 
than the maximum one, J — Jw < Jdw (00), then the disk 
will be completely ablated by the radiation at the radius 
^^out given (from Eq. ^ for = 0) 



(28) 



In this case the outer disk radius i?out is equal to the radius 
above which all material is carried away by the disk wind. 
The corresponding disk wind mass loss rate Mdw(^out) is 
then given by Eqs. (PO]) . Note however that the formu- 
lae discussed in Sect.[5]are strictly valid only in the optically 
thick part of the disk (see Sect. lA.ip . 

If the net angular momentum loss rate J — Jw is larger 
than the maximum one, J—Jw > Jdw (00), then the disk will 
be only partly ablated by the radiation. The net angular 
momentum loss J — Jw is in this case the sum of the parts 
carried by the disk and disk wind. 



j - Jw = jdw(oo) + ja{Md), 



(29) 



where Ja{Md) is given by Eq. p3|) . Here one can assume 
a conservative estimate of the isothermal disk with p — 0. 
From Eq. the mass-loss rate carried away purely by the 
disk Md can be calculated, giving the total required mass- 
loss rate as a sum of parts carried finally by the stellar wind 
(Mw), disk wind (Mdw), and purely by the disk (Md) as 



M = Mw + Mdw(oo) + Md. 



(30) 



The calculation of the functions Jdw{r) and Mdw('') re- 
quires the knowledge of the line force parameters Q and a. 
As the NLTE calculation of these parameters for the disk 
wind environment are not available, one can use their val- 
ues derived for line driven winds for solar metallicity, i.e. , 
Q fti 2000, a nd a « 0.6 (iGavlevI [iQQSl iPuls et al.l ' l2Ciool 
iKrtickal |2006[ ) . For the metallicities other than the solar 
one the scaling Q ~ ^ can be used (here Z is the mass 
fraction of heavier elements) , which is in a good agreement 
with the resu lts of NLTE wind models (jvink et al.ll200lL 
IKrtickal [2nn6h . 



7. Other processes that may influence the outer 
disk radius 

In addition to the radiative force, there may be other pro- 
cesses that may influence the outer disk radius and conse- 
quently determine the required mass-loss rate for a given 
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angular momentum loss rate. For example, in binaries the 
outer disk edge may be naturally truncated due to the pres- 
ence of the companion. However, the most uncertain part 
of the proposed model is connected with the mechanism of 
the viscous transport, which may also influence the outer 
disk radius. 

7.1. Loss of viscous coupling 

The magnetorotational instability (jBalbus fc Hawlevll 19911 ) 
is a promising mechanism to explain the source of anoma- 
lous viscosity in accretion disks. As the dynamics of accre- 
tion and decretion disks is similar, it is likely to be impor- 
tant also for the angular momentum transfer in decretion 
disk. However the nu merical s imulations of magnetorota- 
tional instability (e.g. IStone et al. 1996 , Hawlev & Kroli^ 
concentrate on the inner parts of the disk, whereas 
the evolution close to the sonic point is, to our knowledge, 
not very well studied. The stability con dition of the posi- 
tive d erivative of the angular frequency (jBalbus fc HawlevI 
[T99I d^Vd r > is fulfilled even in the supersonic wind 
region. The fact that the ratio of the viscous timescale 

''vise ~ {^aVl (H/r)^^ to the growth timescale of the mag- 
netorotational instability tmri « 1/^2 decreases with radius 
as TvIsc/tmri ~ {vk/o.) /a (jHavasaki fc Okazakill200d ) in- 
dicates that in the outer parts of the disk where the az- 
imuthal velocity is lower than the thermal speed the mag- 
netorotational instability would not be effective. As this 
happens at supersonic velocities, this again supports our 
conclusion that Eq. (jlSp indeed gives the upper limit for 
the angular momentum loss. 

Moreover, the ratio of the particle kinetic energy to 
the absolute value of its gravitational potential energy is 
roughly equal to 1/4 at the critical point. Consequently, 
for radius few times larger than the critical one the disk 
material may freely escape the star and the viscous sup- 
port is no longer needed. 

The loss of the viscous coupling may occur even before 
the radial disk expansion becomes supersonic. In such a 
case for a — >■ from Eq. follows that the disk starts 
to be momentum conserving and the location of the point 
where this occurs sets the outer disk radius i?out- 

Note also that the disk equations were derived assuming 
that the disk is geometrically thin, i.e., H <^r. However, at 
the critical point the ratio of the disk scale height to radius 

H/r « \J~^ is of the order of unity (assuming isothermal 

disk) . Consequently, above the critical point the vertical av- 
eraging used for the obtaining of Eqs. ([T])-® is no longer 
applicable. On the other hand, because the angular mo- 
mentum loss reaches a plateau below this point, this effect 
has not a significant influence on our results. 

7.2. Influence of stars in a close neighbourhood 

For members of binaries or for stars in a very dense star 
cluster the disk can be potentially truncated due to the 
influence of a nearby star. 

The nearby star could disrupt the disk by its gravita- 
tional interaction with the disk. In this case the outer disk 
radius is that at which the gravitational field of the nearby 
star starts to dominate, i.e., the Roche lobe radius in the 
case of binaries. 



If the disrupting star is luminous one, then it may dis- 
rupt the disk via the radiative force. This case is analogous 
to the case of the radiative ablation due to the central star. 
Consequently, we conclude that this effect would be impor- 
tant only if the nearby star is located within a few radii 
from the central star. 

Finally, the nearby star may heat the disk material in- 
creasing the local sound speed, and consequently decreasing 
the critical radius above which the disk material may leave 
the star. 

Taking all discussed disruption mechanisms together, 
we conclude that in the case of the nearby star with a sim- 
ilar spectral type the disruption is effective only if the dis- 
rupting star is at the distance lower or comparable to the 
critical radius. If the nearby companion is able to disrupt 
the disk, the angular momentum loss becomes less efficient, 
and the star has to lose a larger amount of mass to keep 
the rotation velocity below the critical one. Consequently, 
we expect larger disk mass-loss in close binaries and in very 
dense star clusters. 



7.3. The disk build-up and its angular momentum 

In the analysis presented here we used an assumption of 
constant required angular momentum loss rate, which en- 
abled us to use stationary equations. This assumption is 
reasonable in most phases of the stellar evolution, as the 
evolutionary timescale is much longer than the typical 
timescale of the disk build-up, wh ich is of the order of years 
()Okazakill200l.lJones et al.ll2008D . This also means that the 
transitional processes that occur when the star reaches or 
leaves the critical limit are more complicated than studied 
here. 

In the course of the stellar evolution, when the surface 
rotational velocity reaches the critical limit, in a first time 
the disk appears because it is feeded by the mechanical 
mass loss. The disc grows and part of it is ablated and part 
is transported away via viscous coupling until an equilib- 
rium between the required angular momentum loss rate and 
mass-loss rate is achieved. During this process the disk own 
angular momentum could be of some importance. 

On the other hand, when the star leaves the critical 
limit, an inner part of the remaining disk is accreted on 
the star while other par ts are expelled into the interstellar 
medium (|Okazakill200l . 

7.4. Implication for stars with disk 

The processes discussed here might be relevant also for 
other stars with disks. For example, the disk radiative ab- 
lation might be one of the reasons why the Be phenomenon 
is typical for B stars only, whereas for more luminous O 
stars any disk could be destroyed by the radiative force. 

Similar effects should also be present in accretion disks 
during star formation. In more luminous stars the radia- 
tive ablation could contr i bute to the disk photo evaporation 
fe.g-. lAdams et al.|[200l [Alexander et al.ll2006[ ) in dispers- 
ing of the disk. Moreover, a similar process of the angular 
momentum transfer is present also in the accretion disks of 
these stars, consequently influencing the distribution of the 
rotational speeds on the ZAMS. 
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7.5. Implications for first stars 

Mechanical mass loss through a decretion disk can be a 
ubiquitous phenomenon especially for Pop III or very metal 
poor stars. Indeed as shown by Ekstrom ct al. (2008a) pure 
hydrogen- helium Pop III stars with masses above 60 Mq, 
beginning their evolution on the ZAMS with a surface ve- 
locity around 70% of the critical angular velocity, will reach 
the critical velocity during the MS phase. This arises be- 
cause of two effects: first angular momentum is transported 
from the inner regions to the surface during the MS phase; 
second, the angular momentum accumulates at the sur- 
face since it is not removed by stellar winds. Note that 
in the absence of metals hydrogen and helium are unable 
to drive a line-driven wi nd being nearly completely ionized 
(iKrticka fc KubatJl200fil^ . 

As hydrogen-helium first stars are unable to launch a 
line-driven wind, we expect the radiative ablation to be 
inefficient close to the star. On the other hand, at larger 
distances a nonnegligible fraction of hydrogen could become 
neutral, enabling the possibility of disk radiative ablation. 

The disk wind mass-loss rate in such c_ase could be de- 
scribed as a flow with a very low value of Q (corresponding 
likely just to Lya line force). A rough estimate of the disk 
wind mass-loss rate in this case could be obtained insert- 
ing instead of Mcak the single line mass-loss rate estimate 
M K. L/c^ (jLucv fc Solomonlll970[ ) in Eq. Anyway, in 
most cases such flow would be likely inefficient, especially 
because the disk wind mass-loss rate originates close to the 
star (see Fig. IA.2I) . Consequently, the relation between the 
mass-loss rate required for a given angular momentum loss 
rate would be given by the wind-free condition Eq. (1151) . 

7.6. Future work 

The most uncertain ingredients of a proposed model are the 
viscous coupling, the disk temperature distribution and the 
radiative ablation. To include these processes we applied 
the same description used in the accretion disk theory and 
the theory of radiatively driven winds of hot stars. This may 
not be completely adequate for the description of decretion 
disk especially at large distances from the star studied here. 
Consequently, future work should address these problems. 



8. Conclusions 

We examine the mechanism of the mass and angular mo- 
mentum loss via decretion disks associated with near- 
critical rotation. The disk mass loss is set by the angular 
momentum needed to keep the stellar rotation at or below 
the critical rate. We study the potentially important role 
of viscous coupling in outward angular momentum trans- 
port in the decretion disk, emphasizing that the speciflc 
angular momentum at the outer edge of the disk can be 
much larger than at the stellar surface. For a given stellar 
interior angular momentum excess, the mass loss required 
from a decretion disk can be significantly less than invoked 
in previous models assuming a direct, near-surface release. 

The efficiency of the angular momentum loss via disk 
depends on the radius at which the viscous coupling ceases 
the transport the angular momentum to the outflowing ma- 
terial. When the radiative force is negligible, we argue that 
this likely happens close to the disk sonic (critical) point 



setting the most efficient angular momentum loss. In the op- 
posite case, when the radiative force is nonnegligible, there 
is not a single point beyond which the viscous coupling dis- 
appears. The disk is continuously ablated below the sonic 
point, and the ablated material ceases to be viscously cou- 
pled, decreasing the efficiency of angular momentum loss. 

We describe the method to include these processes into 
evolutionary calculations. The procedure provided enables 
to calculate the mass-loss rate necessary for a required an- 
gular momentum loss just from the stellar and line force 
parameters. We can distinguish three different physical cir- 
cumstances: 

case A: When the disk wind is able to remove the whole ex- 
cess of angular momentum (the disk is completely 
ablated by the wind, see Eq. ([25]) ) then the outer 
disk radius is given by Eq. (I28|) . and the required 
mass loss is given by Eq. ([20)1 . The limiting case 
-Rout ~ Rcci would then correspond to the near 
surface release of the matter without any disk. 
Note that in a rare case when the analysis leads 
to i?out > ^crit the radius i?crit should be used 
as the outer disk radius (case B). The expressions 
presented in the paper are given in the hypothesis 
of an optically thick disk and should be appropri- 
ately modified for optically thin disks. 

case B: If the raditiave force is not able to remove sufficient 
angular momentum (the disk is not completely ab- 
lated) then part of the excess angular momentum 
must be carried away by the disk (Eq. ([^ V In 
this case Eqs. (^5]) . ([201) can be used to estimate 
the mass-loss rate. The outer disk edge could be 
identified with the critical point. 

case C: If the effects of the radiative force are negligible, 
then the whole excess of angular momentum is car- 
ried away by the disk and the the outer disk edge 
is approximately given by -Rcrit and the required 
mass- loss rate could be derived from Eq. ((T5|) . 

Finally, we note that, in absence of strong magnetic 
field, many of the features discussed here may also be ap- 
plicable to the case of star-formation accretion disks. 
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Appendix A: Disk wind mass-loss rate 

A.l. Disk optical depth 

In the case when the disk is optically thick in continuum, 
the disk outflow may be driven not only by the radiation 
from the stellar surfac e, but also by the st ellar radiation re- 
processed by the disk (jGavlev et al.lll999l ) . To estimate the 
optical depth of the disk, let us assume hydrogen and he- 
lium to be ionized in the disk. In this case a significant part 
of the disk optical depth originates due to the light scatter- 
ing on free electrons (for wavelengths lower than that cor- 
responding to the Balmer or Lyman jump also bound-free 
transitions may contribute). The transverse optical depth 
is then roughly given by r = J Kepdz = Ke^, where Kc is 
the Thomson scattering cross-section per unit of mass. The 



disk is optically thick in the vertical direction (t > 1) if the 
mass-loss rate is larger than 



M > 



2-KrVr 



10"^^ Mq year" 



IRr. 



Ums-1 



(A.l) 

For a given mass-loss rate the disk is optically thick close 
to the star, while becoming optically thin at larger dis- 
tances. For example, for a typical disk mass-loss rate re- 
quir ed by the evolutionar y calculations 10~^ M0/year~^ 
(e.g. lEkstrom et al.l[2008b[ ) the disk is optically thick even 
at large distances from the star r w lO^R© for subsonic 
radial velocities. Consequently, in realistic situations the 
disk is likely to be optically thick, at least close to the 
star, resembling the " pseudophotosphere" of Be stars (e.g. 
iKoubskv etaDllQQTl) . 

Contrary to very dense hot star winds (where the radia- 
tive flux comes from regions below the photosphere), here 
we expect that the wind from the optically thick disk starts 
to accelerate above the point where the disk optical depth is 
unity. Numerical results show that the height of this point 
is comparable to the disk scale height H for moderate disk 
mass-loss rates M < 10~^ M0/year~^. Consequently, we 
shall neglect the disk geometrical height in in our analyze 
here and assume that the disk wind originates from the 
equatorial plane. 



A. 2. Disk wind equations 

The outflow from the optically thick disk irradiated by the 
central star can be understood within the framework of th e 
wind driven by external irradiation (|Gavlev et al.l |1999[ ) . 
We study the disk outflow in noninertial frame corotat- 
ing with the disk. We use the Cartesian coordinates with z 
axis perpendicular to the disk (see Fig. lA.ll) . The disk wind 
originates in the disk plane z — 0. We assume purely ver- 
tical flow with velocity Vz{z) and we neglect a potentially 
important part of the radiative force due to the Keplerian 
velocity gradient (jGavlev et al.ll200lh . 
The stationary continuity equation 

V (pv) = 0, (A.2) 

takes within our assumptions the form of 

d 

— {pvz) = 0, or m = pvz — const., (A. 3) 

oz 

where m is the disk wind mass-loss rate per unit of disk 
surface. 

The radiative force of the ensemble of lines in the 
Sobolev approximation is then ((Rvbicki fc Hummeil I1978L 
ICranmer fc Owockilll995l . lGavlevlll995'l) 



-2a 



Srad 



I- a 



/(n) 



/nV (nv) 
V P 



ndrj, 



where Q and a are line force parameters. Ignoring the in- 
coming beam, and simply assuming all the locally normal 
incident radiation from one hemisphere is directly reflected 
upward in vertical beam normal to the disk, the intensity 
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Fig. A.l. The coordinate system for the calculation of the radiative force 



is given by 



COS0 / d/iyl — /i^/* 



1)^^ / ^1 - ^'dM, (A.5) 



where J* = (r/R) F/tt is the emergent intensity from the 
stellar photosphere, F is the radiative flux at radius r, R 
is the stellar radius, ji = cos9, cj) are spherical coordinates 
with origin at the stellar centre, /i* — a/1 — R^ /r"^, and /i, 
and are the direction cosine and azimuthal angle mea- 
sured from the disk plane (see Fig. IA.l| ). The z-component 
of the radiative force in this case is (Eq. (jA.4[) . lG"avlev et al] 
1199911 



5rad — C 



dz 



(A.6) 



where 
/. = 
and 



7^ 



arccos(/x* 



C = 



1 - a 



KqQ 



l-a 



F r 

TT R 

(A.7) 
(A.8) 



The z-component of the momentum equation including 
the gravity term and neglecting the gas pressure term is 



■ dz 



C 



dvz 
dz 



GMz 



(r2 + z2) 



3/2 ■ 



(A.9) 



The vertical momentum equation Eq. (jA.9|) can be solved 
using the transformations 



where the subscript c denotes the critical point values, from 
which using Eq. (|A.11|) we derive 



Cc 



1 - « (1 + 



^3/2- 



(A.13) 



The location of the critical point above the disk plane can 
be derived from the regularity condition (CAK), which 
yields that the critical point occurs at the point of max- 
imum of z component of the gravity acceleration at a given 
streamline, 

Hence, the point of the maximum acceleration a.cts as the 
throat of the nozzle flow (jFeldmeier fc Shlosmanl[l999l ). 

The total disk wind mass-loss rate is then given by an 
integral of the mass-loss rate per unit of the disk surface m 
between the equatorial radius i?oq and the outer disk radius 

^out 

/■-Rout 

Mdw(-Rout) = 2 X 27r / mrdr, (A.15) 
where from Eqs. (TOI) . (fOOj) . (f02t 



1 / KcQ 



JzO_Y (GM\ 
l-a) \ J 



(A.16) 



The factor of 2 in Eq. (|A.15I) comes from the fact that the 
wind originates from both sides of the disk. Consequently, 
the total disk wind mass-loss rate is 

Mdw(i?out) = ^1 {TQ)-^P^ (^) . (A.17) 
where the Eddington parameter F = Ko-^/ (47rGAi^c), and 



yielding 



c 

K 



2GM 

z 



w' = Kw" 



c 



(1 + C^) 



2^,3/2^ 



(A.lOa) 
(A.lOb) 

(A.lOc) 
(All) 



where the prime denotes the derivative with respect to (. 
This equation has a critical point 



aKcW^ 



fa-l 



0, 



(A12) 



PliXout) = 



l-a 



3/2 



F J 



and F is the flux at radius r. Comparing with the CAK 
mass-loss rate estimate Eq. (1161) we have 



(A.19) 



Mdw(-Rout) = ^^1(^1 McAK. 



Assuming the disk wind is not viscously coupled to the 
disk, the total angular momentum loss via the disk wind is 



Jdw(i?out) = 2 X 27r / r^v^mdr, (A.20) 
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0.08 



a = 0.6 



the value of P\i2 = 0.08 indicates that the angular momen- 
tum loss from the disk wind is also lower than the angular 
momentum loss due to the stellar wind. 



Q, 



0.06 



a = 0.6 



0.02 



0.04 








a = 0.4 



a = 0.4 



10 



100 



x= r/R 



Fig. A. 2. The behaviour of Pi{x) (solid line), and Pl{x) 
(dashed line, Eq. (jA.181) ') for different values of CAK pa- 
rameter a. 

where m is given by Eq. (jA.16p . As the disk wind originates 
mainly from the regions close to the star (with r/R < 10, 
see Fig. lA.2|) . where the azimuthal velocity is roughly equal 
to the Keplerian one (see Fig.[T]), we can assume « VK{r) 
in Eq. (jA.20l) and consequently 



Again, using the CAK mass-loss rate estimate Eq. (IT5|) the 
disk wind angular momentum loss 



is by a factor of Pi {Rout/ R) larger than the angular mo- 
mentum loss due to the CAK wind launched from equator 
of hypothetical critically rotating spherical star with radius 



Expansion of Eq. (jA.7|) about R/r = shows that far 
away from the star, fz/F ~ 2i?/(37rr). The comparison 
between a precise formula Eq. (|A.7p and its asymptotic 
form shows apart from a small region 1 < x < 3/2, the 
agreement is actually quite good. If we use this asymptotic 
form over the full range from x = 3/2 to Xout in Eq. (jA.18|) 
we find (using Eqs. (|XT3|) . (|XT4l) ') 



Our analytical results (s ee Fig. IA.2I) a re in agreement with 
numerical calculations of lProea et al.l 11998) that show the 
disk mass loss is dominated by material arising from the 
inner region of the disk (r < 10 i?). Consequently, if the 
radiative force is not strong enough to disrupt the disk close 
to the star, it is unlikely that it would be able to do so in 
the outer parts of the disk. 

The modern hot star wind models give slightly lower es- 
timate of the mass-loss rate than the CAK formula Eq. (fT6|) 
due to inclusion of finite disk correction. However, be- 
cause formula Eg (IA.23I) for a typical value of a « 0.6 
(|Puls et all I2OOOL iKrticka 2Q0i) gives for Xout ^ 00 the 
value of Pi ~ 0.04, the disk wind mass-loss rate is signifi- 
cantly lower than the stellar wind mass-loss rate. Similarly, 




(A.21) 




(A.22) 



R. 




